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Abstract 

We study the theta dependence of the glueball spectrum in a strongly 
coupled cousin of large N gluodynamics defined via the AdS / CFT correspon- 
dence. By explicitly diagonalizing the lOd gravity equations in the presence 
of the RR 3-form and 1-form fluxes we found a mixing pattern for the lowest- 
spin lightest glueballs. The mixing between the scalar and pseudoscalar states 
is not suppressed, suggesting that the CP-odd effects persist in the large N 
theory. As a consequence, the lightest mass eigenstate ceases to be a parity 
eigenstate. We found the former as a linear combination of a scalar and pseu- 
doscalar glueballs. On the other hand, the mass eigenvalues in a theory with 
and without the theta term remain equal in the large N limit. 



1 Introduction and summary 



QCD in the limit of a large number of colors, oo, is expected to contain a great 
deal of information on confinement ^j. In this limit gluons dominate the partition 
function and confinement should be easier to study in the theory without quarks. 
The subject of the present work is large gluodynamics (i.e., QCD without quarks) 
with a nonzero theta term. Observationally, the theta parameter, if not zero, should 
be truly tiny, 9 < 10~^, and, studies of the theta dependence in the large N theory 
may seem less motivated. However, this is not so because of the following arguments. 
Large gluodynamics with the theta term reveals a rich ground state structure with 
multiple vacua separated by domain walls j2[Ell31in|- This by itself is an extremely 
interesting manifestation of nonperturbative physics in a non-Abelian gauge theory. 
Moreover, certain remnants of the large A^ vacuum structure are expected to be 
present in real QCD too. Examples of this include the heavy nonperturbative states 
that make the domain walls [B], and new properties of axion and hadronic domain 
walls in a theory where the strong CP problem is solved by the axion mechanism 
(for recent summary see Ref. [7j). 

The theta term enters the QCD action suppressed by one power of A^ compared to 
the gluon kinetic term. This may seem to indicate that no theta dependence should 
survive in the large A^ limit. For instance, a dilute instanton gas approximation 
would give rise to the theta dependence that is suppressed as exp(— A^). However, 
there is a body of evidence (lattice and otherwise) suggesting that this is not so. The 
best example is the theta dependence of the vacuum energy [HI I2l El El which is of 
order 0{N^). This dependence can only be attributed to certain infrared effects that 
are not captured by a dilute instanton gas approximation. The related issue that 
has not been explored so far is the theta dependence of the spectrum of glueballs in 
this theory. Masses of glueballs are generated by noperturbative effects as well [TU] . 
Similar effects could in general introduce the theta dependence in the spectrum even 
in the large A^ limit. In particular, we did not find convincing arguments to believe 
that the theta dependence would be completely washed out from the spectrum in 
the large A^ limit. How would one study these issues in more detail? 

The theta term introduces a complex phase in the euclidean formulation of the 
partition function of a theory the action of which is otherwise real and positive 
semidefinite. Because of this lattice studies of these issues seem difficult (as an ex- 
ception, see, e.g., IHl). One way to address the question of the theta dependence 
of glueballs is to use the AdS/CFT correspondence [121 O (HI adapted to non- 
supersymmetric and non-conformal theories (for a review see, ^21). A well-known 
drawback of this approach is that it only allows to calculate gauge theory observables 
for a strongly coupled cousin of large A^ gluodynamics which in addition contains 
states that are not present in QCD. However, in the large A^ limit the effects of 
the additional states on QCD resonance physics should be minimized, and in the 
absence of other methods this seems to be a reasonable starting place. 
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In the main part of the work we study the theta dependence of the glueballs 
by looking at the HA construction of N D4-branes with one compact worldvolume 
dimension JH] and a nonzero bulk RR 1-form flux. The latter gives rise to the theta 
dependence in the dual gauge theory 0. In this setup we manage to diagonalize 
exactly the bulk gravity equations for the fluctuations of a graviton, dilaton and RR 
1-form. The diagonalization of the gravity equations suggest the following remark- 
able pattern for the theta dependence of the gauge theory glueballs. The lowest 
scalar (0^"*") and pseudoscalar (0~^) glueballs mix. There are kinetic as well as 
mass mixing terms both of which are determined by 6. The system of glueballs can 
be diagonalized by shifting the field of the lightest spin- zero state by a field that is 
proportional to the heavier spin-zero glueball multiplied by 6. Interestingly enough, 
the above diagonalization can only be achieved by shifting the field of the lightest 
spin-zero state, which in gluodynamics is necessarily the 0"'""'" glueball The 
heavier state remains intact. This property, as we will discuss in Section 4, is vital 
for the applicability of the method of the calculation. After the diagonalization, the 
lowest mass eigenstate ceases to be a parity eigenstate. The mixing between the 
parity eigenstates is a leading effect in large N (i.e., it is not suppressed by powers 
of N). However, the masses of the physical (diagonal) states in a theory with and 
without the theta term are the same. Hence, the only effect of the theta term in 
the leading large limit is that the lowest mass spin-zero state becomes a mixed 
state of O"*"^ and O""*" glueballs. The periodicity of the wavefunction with respect to 
6^6 + 27r (integer) is achieved by choosing appropriate branches of the theory in a 
way that also makes the vacuum energy periodic in 6' [3 0] • 



In order to study the theta dependence of glueballs from the Gauge/Gravity corre- 
spondence we consider the dual supergravity description of pure U{N) gauge field 
theory introduced in (TH]. In this setup one starts with weakly coupled type HA 
superstring theory in the presence of N Z)4-branes. The D4 soliton background 
solution relevant for this discussion is obtained by compactifying on a circle of ra- 
dius M^]^ one of the 4 spacial directions of the 4-brane solution and by taking the 
near horizon limit. The effect of the compactification in the D4-branes worldvolume 
theory is that the fermions acquire masses of order Mkk at tree level due to the anti- 
periodic boundary conditions and so do the scalars due to loops. The worldvolume 
theory at energies below Mkk is then 4-dimensional nonsupersymmetric noncon- 
formal U{N) gauge theory. Thus we proceed to review the DA soliton solution and 
how the effect of a theta term is incorporated. 

The low energy effective action of type HA superstrings in the Einstein frame 
and with vanishing Kalb-Ramond B field reduces to 



2 L>4 Soliton with RR 1-form flux 



1 



/ 




(1) 



(27r)7/f 



(4) , 
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where Ig and gs are the string length and couphng constant and F(2) and F(4) are 
the field strengths of the RR 1 and 3-forms respectively. From these, the equations 
of motion for the metric, dilaton and 1-form follow: 



UfmustF^^'^ - IgMNF^^ , (2) 



48 

^;^^5m (v^/^"M = r''^'Fm + j-s-'^n)^ (3) 

5a/ (v^e^/'i^*'^) =0. (4) 

The so called DA soliton is given by the following Einstein frame metric, dilaton and 
constant 4-form: 



(^) 



= (^) ^ ^(4) - 9;'e, , 

R' = T^gMl , f{U) = 1 - {UKKlUf , (6) 

where x^, fi = 1 . . .4 are the euclidean coordinates of the noncompact directions 
and r denotes the coordinate along the compactified circle on the DA worldvolume. 
dO| and are the line element and volume form on a unit 4— sphere and U is 
the coordinate on the radial direction transverse to the DA branes. In order to 
avoid a conical singularity at U = Ukk, t is identified with period A = 27rM^^ = 

A7TR^/y3Ul/K- 

The supergravity description is valid in the regime of small curvatures (in string 
length units) and string coupling, namely 

< 1 , gse"^ < 1 . (7) 

The latter implies a maximum value Umax of the radial coordinate that should 

therefore be much larger than the parameter Ukk P^- This regime corresponds 

to the 't Hooft limit of the four dimensional theory gyM 

fixed ^ 1, where the identification used gyM ~ 9sIsMkk can be read off from 

below. 

The theta dependence of this background is obtained by turning on the RR 
1— form C(i) in the r direction jH]: 

- -^Tp ' 

where dc = 6 + 2mT with integer n and the normalization has been chosen such 
that / drdUFjjT- = 9c. Indeed, it is not difficult to see that (jH)) solves (jU if the 
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background geometry is not altered by the presence of this new field. This absence 
of backreaction on the geometry is due to the fact that in the large limit the 
contribution from F(2) to © and (jH)) is subleading because the other terms contain 
extra powers of g^^ ~ A^. 



3 Supergravity modes 

In order to find the glueball spectra, we proceed to study the equations of motion 
to linear order in the background of the previous section. In particular we are 
interested in the scalar modes that act as sources for the scalar and pseudoscalar 
operators of the boundary theory. For the case of vanishing theta this modes have 
been identified in [THj and [201 ■ In this section we will find the particular combination 
that diagonalizes the supergravity equations of motion when theta is turned on. 
For the metric and dilaton perturbations we consider the following linearization 

gMN = 9mn + huN , = 0^ + (50 , (9) 

where gf,jj^ and 0^ stand for the background values and the diagonal components 
of the fiuctuations are parametrized as follows: 

hMM = gl,i,v''H{U)e'^-^' , (10) 
Ji^^^^H{U) . (11) 



Respecting the 5*0(4) and 5*0(5) symmetries of the . . .x^ and 5^ directions there 
are tree scalar fiuctuations T, L and 5 (in the notation of 21 ) that correspond 
to the dilaton, 4— sphere volume fiuctuation and the exotic polarization [20] corre- 
spondingly. To simplify the expressions we choose fc^ = S^'^k^ via an 50(4) rotation. 
These are given by: 

T: v"" = i(l,-|-|-|l,l,l,l,l,l), 00 = ^, (12) 
L: v"" = (-1,-1,-1,-1,-1,-1,1,1,1,1), 00 = -^, (13) 

The exotic polarization 5* has also (in the gauge chosen) off diagonal components 

hu.. = = -'Y^^^^HiUy'^^' . (15) 

For each of these polarizations, the function H = T, L, S respectively satisfies: 

U{U^ - 1)T" + {AU^ - 1)T' - k'^UT = , (16) 
U{U^ - 1)L" + {AU^ - 1)L' - {k^U - 18U)L = , (17) 

U{U' - 1)^" + (4f/^ - 1)^' -u(^e- j^^^^j S = 0, (18) 
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where, for convenience, we have rescaled the coordinates to form dimensionless quan- 
tities, namely 

and dropped the tildes. 

Let us now turn our attention to the fluctuations of Cr- We linearize the equation 
of motion for the RR 1-form with fluctuations given by 

a = + X(t^)e'''^"' , (20) 

where is the background value. 

The key point is that the resulting equation can be decoupled from the metric 
and dilaton by choosing 

X = ^^+^/(f/)(T + 5), (21) 

with V satisfying the following equation: 

U{U^ - 1)V" + 4(f/3 - 1)V' - k'^UV = , (22) 

the equation for fluctuations of Cr in the absence of theta term. Thus, we were able 
to diagonalize the supergravity dual of SU{N) gauge theory in the presence of a CP 
violating term. This has important consequences that we discuss in the following 
section. 



4 Couplings to boundary theory 

We infer the coupling of the supergravity modes to gauge invariant operators of the 
four dimensional field theory by considering the Born-Infeld action describing the 
low energy worldvolume excitations of the D4 branes. The operators of interest 
for us are f^^f'^'^ and f^uf'^'^ for Yang-Mills field-strength f^^ and its dual J^'^. 
Neglecting worldvolume scalars it reads (in string frame): 

Ibi = %Tt J d'x (^e" ^det (w + 27ra'/„„) + ^2(27ra')'e™*C^/„./.t) , (23) 

x"^ = r, . . .x^ and T4 = {2TT)^'^g^^lj^ . Upon compactification and expansion to 
second order the couplings to the gauge invariant operators are obtained [IH]. We 
give the expression in Einstein frame: 

Ibi = jU2na'yTT J d^x {e-^^"^r^deig,,f + \iCrff) ■ (24) 
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By neglecting the metric fluctuations proportional to A;^ that give no contri- 
bution when contracted with the conserved energy momentum tensor of the four 
dimensional field theory we get the following result: 

Ibi = IEi+ f d'x UO, + xOA , (25) 



where I^j accounts for the contribution from the background metric and dilaton 
fields and where we have defined the scalar and pseudoscalar glueball operators O4 
and O4 and the corresponding couplings ip and x follows: 



3zA 



^4 = ; f,.r , (26) 



ip = T + S, (28) 
X = Iv + ^f^. (29) 



Note that the S'4— volume scalar fluctuation L decouples. 

The AdS/CFT prescription determines the generating functional for the 4- 
dimensional field theory in terms of the low energy effective string theory partition 
function for on-shell fields 

^Jd^x^oO^ _ g-/sGM ^ (30) 

in which the r.h.s. containing the supergravity action Isg is used as an approximation 
to the string theory partition function and ipo are the boundary values of the on-shell 
normalizable supergravity modes (p acting as sources for the field theory operator 
O. Therefore, to obtain the glueball spectrum we can consider the variation of the 
r.h.s. of (j^UI) with respect to boundary sources to obtain two point functions. For 
example, for pseudoscalar glueballs we would consider the following variation 
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-isG[v,i>] 



Vo=0,i/'o=0 



{d,{x)d,{y)) . (31) 



On the other hand, for scalar glueballs we find a mixing. Not only both scalars T 
and S in the combination ip act as a source for O4 but also, when 6c is nonvanishing, 
ip sources O^. The corresponding two point function is obtained from the following 
variation: 

^ ^ 5 \ / 5 0. 6 \ -isaivM\ = (^o,ix)0,iy)) , (32) 



\S^{jo{x) TT SVo{x) J \6ilJo{y) n 6Vo{y) j 

since f{U) 1 at the boundary^. 

^Eq. H32|l makes the values 9 — ±7r special. This indeed should be so since one expects 
spontaneous CP violation to take place for 9 — zLtt |21 El ■ In our approach the form of Eq. H32|l is 
an artifact of the normalizations that we choose to use for the bulk fields. 
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The above results can be summarized as follows. The field ip sources a linear 
combination of the operators O4 and O4, while the field V sources only O4. Hence, 
the diagonalization of the bulk gravity equations in terms of ip and V suggests the 
following mixing pattern for spin-zero glueballs: The lowest state is a linear super- 
position of a former scalar and pseudoscalar states while the heavier state is intact 
(it is just a pseudoscalar state). The mass of the pseudoscalar state, determined by 
Eq. (j22I), is the same as in the 9 = theory (determined in Refs. [Tnil2I|)- So is the 
mass of the mixed lightest state. Does this pattern have any special meaning? It 
seems it does. To see this consider the following two two-point correlation functions 

cix) = {o,{x)o,m - (O4(0))^ 

and 

Cix) = {d,ix)d,io)) - {omf- 

These correlators can be saturated by the corresponding orthonormal set of physical 
intermediate states. Then, the leading behavior of C{x) and C{x) at large euchdean 
is determined by the corresponding lightest states. In a theory without the theta 
term those states are scalar and pseudoscalar glueballs respectively. However, once 
the theta term is switched on the scalars and pseudoscalars can mix. In general this 
mixing could be arbitrary. If so, then the former scalar state would also contribute to 
the expression for C{x). However, because the scalar is lighter than the pseudoscalar 
(and after mixing it can only become even lighter), this would mean that C{x) at 
large is completely dominated by the residue of the former scalar state. If this 
were true, we would not be able to determine the properties of the pseudoscalar 
state by calculating C{x). Fortunately, this does not happen. The contribution of 
the physical lightest state (which is the former scalar glueball) is exactly canceled 
in C{x), as suggested by the diagonalization of the gravity equations. 
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